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Abstract 

Multiple p- and s-polarized compound surface plasmon-polariton (SPP) waves at a fixed frequency can 
be guided by a structure consisting of a metal layer sandwiched between a homogeneous isotropic dielectric 
(HID) material and a periodic multilayered isotropic dielectric (PMLID) material. For any thickness of 
the metal layer, at least one compound SPP wave must exist. It possesses the p-polarization state, is 
strongly bound to the metal/HID interface when the metal thickness is large but to both metal/dielectric 
interfaces when the metal thickness is small. When the metal layer vanishes, this compound SPP wave 
transmutes into a Tamm wave. Additional compound SPP waves exist, depending on the thickness of 
the metal layer, the relative permittivity of the HID material, and the period and the composition of 
the PMLID material. Some of these are p polarized, the others being s polarized. All of them differ in 
phase speed, attenuation rate, and field profile, even though all are excitable at the same frequency. The 
multiplicity and the dependence of the number of compound SPP waves on the relative permittivity of 
the HID material when the metal layer is thin could be useful for optical sensing applications. 


1 Introduction 


The electromagnetic fields of a surface-plasmon-polariton (SPP) wave guided by the planar interface of a 
homogenous metal and a homogeneous isotropic dielectric material decay along the direction normal to 
the interface in both materials [^. The SPP wave is the classical counterpart of a train of quantum quasi¬ 
particles formed jointly by plasmons and polaritons bound to the metal/dielectric interface. The SPP wave is 
p-polarized. Only one SPP wave can be guided by the interface at a specific optical frequency. Due to the field 
being bound tightly to the metal/dielectric interface, the SPP wave is strongly sensitive to the constitutive 
parameters of the materials on either side of the interface, leading to optical sensing applications [2j|^. For the 
same reason, researchers have also focused on surface-plasmon-resonance microscopes with subwavelength 
resolution to obtain images of objects with very low contrast, without the use of dyes or markers [^[^. 
Applications of SPP waves are on the horizon in the area of communications as well [^. 

The excitation of a multiplicity of SPP waves at a specific frequency is very appealing therefore. This 
multiplicity is possible if the isotropic dielectric partner is periodically nonhomogeneous in the direction 
normal to the interface, which has been established both theoretically |8j|^ and experimentally [lOpd] . Even 
though all SPP waves are excited at the same frequency, they will differ in one or more of the following 
attributes: polarization state, phase speed, attenuation rate, and field profile. 

In the present paper, we theoretically formulate and analyze a boundary-value problem that involves two 
parallel metal/dielectric interfaces in a bid to coalesce the SPP waves guided separately by each interface into 
compound SPP waves. For simplicity, one of the two dielectric materials is taken to be homogeneous while 
the other dielectric material is periodically multilayered normal to the interfaces. Both dielectric materials 
are isotropic as also is the metal separating them. 

The plan of this paper is as follows. Section presents the formulation of the boundary-value problem, 
wherein one half space is occupied by a homogeneous isotropic dielectric (HID) material, another half space 
is occupied by a periodic multilayered isotropic dielectric (PMLID) material, and the two half spaces are 
separated by a thin metallic layer, as shown in Figure In this structure, compound SPP waves of the p- 
polarization state can not interact with compound SPP waves of the s-polarization state, thereby splitting the 
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Figure 1: Schematic of the boundary-value problem. A metal layer of thickness L separates a half space 
occupied by a HID material and another half space occupied by a PMLID material of period 217. 


boundary-value problem into two autonomous sub-problems. Numerical results showing the compounding 
of the SPP waves guided by the metal/HID and metal/PMLID interfaces are presented in relation to the 
thickness of the metal layer in Section Concluding remarks are presented in Section 

An exp {—iuit) dependence on time t is implicit, with ui denoting the angular frequency and i = -y/—1. 
Furthermore, kg = LOy/Soiio, Cg = w/Aq, rjg = ^[XgjEg and Ac = 27r/A:o, respectively, represent the wavenumber, 
phase speed, characteristic impedance, and wavelength in free space, with [ig as the permeability and £„ as 
the permittivity of free space. Vectors are in boldface, column vectors are in boldface and enclosed within 
square brackets, and matrices are twice underlined and square bracketed. The asterisk denotes the complex 
conjugate. Cartesian unit vectors are denoted by u^,, u^, and Uj,. 

2 Theoretical framework 

The schematic of the boundary-value problem is shown in Figure The half space z < —L is occupied by a 
HID material with real-valued relative permittivity > 0. A metal layer of thickness L and complex-valued 
relative permittivity Em separates the HID material from the PMLID material which occupies the half space 
z > 0. The unit cell of the PMLID material consists of N layers with real-valued relative permittivity Ej > 0 
and thickness dj (j = 1,..., TV), giving the period 2V, = di + d 2 + ■■■ + djsf. 

Without loss of generality, we consider the direction of propagation of the compound SPP wave to be 
parallel to the x axis. Due to the isotropy of all materials, we can consider p- and s-polarized compound 
SPP waves guided by the structure in Figure separately from each other. 
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2.1 p-polarized compound SPP waves 

The electric and magnetic field phasors of a p-polarized compound SPP wave can be represented everywhere 
as 

E(r) = [ex{z)ux + e:,{z)U:,] exp{iqx) 

H(r) = hy{z)uy exp{iqx) 

where q is the complex-valued wavenumber. The functions ej;(z), ez{z), and hy{z) are determined piecewise 
as follows. 

In the layer of the PMLID material, the Faraday and the Ampere-Maxwell equations yield 
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One of the two eigenvalues represents a field decaying as z —)■ oo, the other a field decaying as 
—oo. Since the fields of the compound SPP wave must decay in the half space occupied by the PMLID 
material [9, 13 
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provided that the eigenvalues are labeled so that Re ^In < 0 and, ipso facto, Re ^In > 0. The 
coefficient Cp is as yet unknown. 

In the half space z < —L occupied by the HID material, field representation is of the textbook variety; 
thus, with bp as an unknown coefficient we have 


E(r) = b, 
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H(r) = —bpUy — exp(zga;) exp [iad{z -\- L)] 
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where Ud = +v^j = k^Sd, and Im (ad) < 0. Accordingly, 
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Inside the metal slab, the Faraday and the Ampere-Maxwell equations yield 
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This differential equation yields 

f(p) (0") = exp (-L+) 
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Standard electromagnetic boundary conditions mandate the equalities (0 )] = Iflp) (O'*')] and 
[f(p) (—L“)] = (—L+)]. After enforcing these equalities and substituting Eqs. ([8| and (10) in Eq. (14), 

we get 
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which can be recast as 
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The dispersion equation for p-polarized compound SPP waves then is as follows: 


det 




= 0 . 


kgUd 

nd 

(15) 

Vo . 

(16) 

is as follows: 
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For any solution q of Eq. (0, the ratio bp/cp can be found using Eq. ( |I^ . 

2.2 s-polarized compound SPP waves 


The electric and magnetic field phasors of an s-polarized compound SPP wave can be represented everywhere 
as 

E(r) = ey{z)uy exp{iqx) 

H(r) = [h,;{z)u^ + h^{z)u^] exp{iqx) 


, z G (—oo, oo) 


(18) 


where the functions ey{z), h^^z), and hz(z) are determined piecewise as follows. 

In the layer of the PMLID material, the Faraday and the Ampere-Maxwell equations yield 
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In the half space z < —L occupied by the HID material, with bs as an unknown coefficient we have 
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Satisfaction of standard electromagnetic boundary conditions leads to the relation 
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The dispersion equation for s-polarized compound SPP waves then is as follows: 
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For any solution q of Eq. (351, the ratio bg/cs can be determined using Eq. (34). 
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3 Numerical results 


Equations ( [T7| ) and (35) were solved numerically using Mathematica (Version 10) on a Windows XP laptop 
computer. For all calculations, we fixed Aq = 500 nm. The HID material and the metal were taken to be 


SFll glass {sd = 3.2508) for the results presented in Secs. 3.1 and 3.2 but Sd G [1,4] was kept variable for the 
results presented in Sec. 3.3 The metal was taken to be silver {e„i = —7.4726 + zO.8). The PMLID material 


was taken to have a unit cell comprising N = 9 lossless dielectric layers of silicon oxide and silicon nitride in 
different ratios [l3). The sequence of the relative permittivities of layers in the unit cell of the PMLID material 
is shown in Table [H While dj = 75 nm Vj S [1, V], the thickness L of the silver layer was kept variable from 
1 to 200 nm. At the chosen wavelength, the skin depth in silver is Sm = 1/Im {kpy/e^) = 29.07 nm. 

In this section, we present the solutions q = qko of the dispersion equations (17) and (35). As our 
computer code had significant errors in the computation of the transfer matrix of the PMLID unit 

cell for Re{q) > 3.0, the data presented in this section are restricted to Re{q) < 3.0. We also present the 
spatial prohles of the Cartesian components of the time-averaged Poynting vector 


P (cc, z) = (l/2)Re [E(a;, z) x H*(a:, z)] 


(36) 


for representative solutions. These profiles were computed by setting either bp = 1 or bg = 1, as appropriate. 

Let us note here that chosen structure can guide surface waves even in the absence of the metal layer, such 
waves being called Tamm waves 14 -18 . Tamm waves have not only been experimentally observed 
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but have also been used for optical sensing applications 21 22 . Our theoretical framework accommodates 
Tamm waves by the simple expedient of setting L = 0. The solutions of Eqs. 0 and ( |35[ ) for Tamm waves 
guided by the HID/PMLID interface when L = 0 are also provided in this section. 


Table 1: Relative permittivity Sj at Aq = 500 nm of the dielectric layer in the unit cell of the 
PMLID material 


j 

1 

2 

3 

4 

5 

6 

7 

8 

9 


4.0720 

3.3787 

3.2176 

3.0628 

2.7936 

2.6488 

2.4759 

2.3388 

2.1986 
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3.1 p-polarized compound SPP waves 

The normalized wavenumbers q calculated for p-polarized compound SPP waves are reported in Table for 
different thicknesses L G [1,200] nm of the metal layer. The table also lists the values of q for (a) the sole 
p-polarized SPP wave guided by the metal/HID interface by itself and (b) the four p-polarized SPP waves 
guided by the metal/PMLID interface by itself. 

Table shows that as many as five different p-polarized compound SPP waves can be guided by the 
HID/metal/PMLID structure for L G [1,200] nm. These waves are organized in five sets labeled pi to p^. 
Compound SPP waves in the sets Pz,Pa, and ps do not exist if the thickness L of the metal layer is < 40.9 nm, 
44.5 nm, and 46.8 nm, respectively. Due to significant computational errors, solutions in the set pi could 
not be found for L < 75 nm. 
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Table 2: Values of q computed for p-polarized compound SPP waves in relation to L. These values are organized in 5 sets 
labeled pi to p^- Solutions obtained for SPP waves guided by either the metal/PMLID interface alone or the metal/HID 
interface alone are also provided, along with the solution for a Tamm wave guided by HID/PMLID interface when L = 0 nm. 
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3.1.1 Superset Va 

The sets pi to ps in Tablej^can be grouped into two supersets. Superset Va comprises the sets pi, pa, p 4 , and 
P 5 . As L increases, the wavenumber of a p-polarized compound SPP wave in any of the four sets in Va tends 
towards the wavenumber of a p-polarized SPP wave that can be guided by the metal/PMLID interface all by 
itself. As L decreases, both the phase speed Vph = Co/Re(g) and the propagation distance = l/AoIin(g) 

decrease in set pi; both Vph and Aprop increase in the sets pa and p 4 ; while Vph is practically constant and 
Aprop increases in the set ps. 

Figure shows the spatial variations of the Cartesian components Px{Q,z) and Pz{0,z) of the time- 
averaged Poynting vector in the plane a; = 0 for representative compound SPP waves belonging to the four 
sets in Va- Whereas Px{0,z) can be discontinuous across bilateral interfaces due to the discontinuity of 
Bz defined in Eq. ([^ across those interfaces, P^(0, z) must be a continuous function of z. These charac¬ 
teristic features are evident in Fig. which also demonstrates that the compound SPP waves are bound 
predominantly to the metal/PMLID interface. The power density is confined mostly to the PMLID material: 
specifically, the first unit cell for compound SPP waves in sets pi and p 3 , and more than one unit cell for 
compound SPP waves in sets p 4 and ps. The affinity to the SPP wave guided by the metal/PMLID interface 
by itself is particularly pronounced when the ratio L/Sm is large, indicative of weak coupling between the 
two metal/dielectric interfaces. 

As L decreases, the compound SPP wave belonging to the set pi becomes also bound to the metal/HID 
interface, thereby indicating a significant coupling between the two metal/dielectric interfaces. However, 
compound SPP waves belonging to pa, P 4 , and ps cease to exist when L/dm ^ 1-4 rather than divert 
significant energy from the PMLID material to the HID material. 


3.1.2 Superset Vb 


Superset Vb comprises only the set p 2 - As L increases, the value of q approaches the value of q for the sole 
SPP wave guided by the metal/HID interface all by itself. However, as L decreases, the value of q approaches 
the value of q for a p-polarized SPP wave guided by the HID/PMLID interface (i.e., L = 0). The set p 2 thus 
exemplifies a transition, and therefore the inherent commonality, between two surface-wave phenomenons 
that otherwise appear to be completely different from each other, and adds to the previously found unity of 
Fano waves and Tamm waves |23|. Furthermore, both Vph and Aprop increase in the set p 2 as L decreases. 


Figure shows the spatial variations of Px{0,z) and Pz{0,z) for representative compound SPP waves 
belonging to P 2 C Vb- Again, while Px{0,z) is discontinuous across bilateral interfaces, Pz(0,z) is a contin¬ 
uous function of z. When the ratio L/Sm is large, the compound SPP wave is predominantly bound to the 
metal/HID interface. As L/Sm decreases, the compound SPP wave also gets bound to the metal/PMLID in¬ 
terface and the energy begins to be transferred from the HID material to the PMLID material. Of course, for 
L — 0 the surface wave is bound to the HID/PMLID interface and is transformed into a Tamm wave HPP)- 


3.2 s-polarized compound SPP waves 

The normalized wavenumbers q calculated for s-polarized compound SPP waves are reported in Table 
for different thicknesses L G [30, 200] nm of the metal layer. The table also lists the values of the three 
s-polarized SPP waves guided by the metal/PMLID interface by itself and the value of the Tamm wave 
guided by the HID/PMLID interface. No s-polarized SPP wave can be guided by the metal/HID interface 
by itself (^[^. 

According to the table, as many as three different compound s-polarized SPP waves can be guided by 
the HID/metal/PMLID structure for L G [30,200] nm. These waves are organized in three sets labeled si 
to S3. Compound s-polarized SPP waves do not exist if the thickness L of the metal layer is < 29.3, 39.3, 
and 43.8 nm for the sets si, S2, and S3, respectively. Therefore it is not surprising that, unlike for compound 
p-polarized SPP waves, the normalized wavenumber q in none of the three sets of solutions approaches the 
value of q for the Tamm wave as L decreases. 


9 






Figure 2: Spatial variations of Poc{0,z) (blue slide lines) and Pz{0,z) (red dashed lines) with respect to z 
for the compound p-polarized SPP waves in Superset Va- (a) L = 75 nm or (b) L = 200 nm in set pi C Va', 
(c) L = 45 nm or (d) L = 200 nm in set ps C Va', (e) L = 45 nm or (f) L = 200 nm in set p 4 C Va', and (g) 
L = 50 nm or (h) L = 200 nm in set p^ C Va- 
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Figure 3: Spatial variations of Px{0,z) (blue slide lines) and Pz{0,z) (red dashed lines) with respect to 0 
for the compound p-polarized SPP waves in Superset Vb- (a) L = 0 nm, (b) F = 1 nm, (c) L = 50 nm, or 
(d) L = 200 nm in set P 2 - 
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Because the metal/HID interface by itself cannot guide s-polarized SPP waves, the wavenumber of the 
compound SPP waves in any of the three sets tends towards the wavenumber of an s-polarized SPP wave 
guided by the metal/PMLID interface by itself as L increases. Finally, as L decreases, the phase speed Vph 
decreases but the propagation distance Ap^-op increases in all three sets Si to S 3 . 

Figure shows the spatial variation of the Cartesian components Px{0,z) and Pz{0,z) of the time- 
averaged Poynting vector in the plane cc = 0 for representative compound SPP waves belonging to the three 
sets of solutions. Unlike for compound p-polarized SPP waves, both components are continuous functions 
of z due to the continuity across interfaces of ej,(z), hx{z) and, in particular, hz{z) by virtue of Eq. ( [2^ . 
The figure shows that the compound s-polarized SPP waves are bound almost totally to the metal/PMLID 
interface resulting in a negligible coupling between the two metal/dielectric interfaces. As a consequence, the 
power density is confined to the PMLID material: specifically, the first unit cell for compound SPP waves 
in the set si, and more than one unit cell for compound SPP waves in the sets S2 and S3. 


Table 3: Values of q computed for s-polarized compound SPP waves in relation to L. These 
values are organized in 3 sets labeled si to S3. Solutions obtained for SPP waves guided by 
the metal/PMLID interface alone are also provided, along with the solution for a Tamm wave 
guided by HID/PMLID interface when L = 0 nm. 


L (nm) 

Si 


S2 


S 3 

30 

1.7375+ f8.13 X 

10"® 

- 


- 

40 

1.7347+ f6.85 x 

10-4 

1.5403-fi6.44 x 

10-4 

- 

50 

1.7335-fi8.91 x 

10-4 

1.5392-ff6.02 X 

10-4 

1.3066-ff3.87 X 10-4 

75 

1.7327-fi9.93 x 

10-4 

1.5384-ff9.37 x 

10-4 

1.3059-ff8.95 X 10-4 

100 

1.7326 + f9.98 x 

10-4 

1.5383-fi9.74 x 

10-4 

1.3057-fi9.62 x 10-4 

150 

1.7326 + f9.98 x 

10-4 

1.5383-fi9.78 x 

10-4 

1.3057-fi9.71 X 10-4 

200 

1.7326 -fi9.98 x 

10-4 

1.5383-fz9.78 x 

10-4 

1.3057-ff9.71 X 10-4 

Qmet/PMLID 

1.7326 + f9.98 x 10"^ 

1.5383-fi9.78 x 10-4 

1.3057-fi9.71 X 10-4 

Qhid/pmlid 



1.8375 




3.3 Effect of the relative permittivity of the HID material 

The data provided in Secs. |3.1| and |3.2| demonstrate that when the thickness L of the metal layer is small 
enough the two dielectric/metal interfaces are significantly coupled each other and compound SPP waves 
propagate bound to both interfaces. Moreover, only one compound SPP wave (in set P 2 ) exists for L smaller 
than the skin depth dm in the metal. But this paucity of compound SPP waves when L < dm not a 
universal feature, as became clear when we fixed L = 25 nm but varied the relative permittivity Sd of the 
HID material from 1 to 4. 

Figures and respectively, provide plots of Re( 9 ) and Ap^op as functions of Ed- A multiplicity of 
compound SPP waves is evident. We have organized these waves in sets labeled pi to ps and si to S3, 
according to the polarization state. Six compound SPP waves, three p polarized and three s polarized, 
propagate for Ed ^ 1.3. These waves begin to disappear one by one ss Ed increases, reducing to a sole 
p-polarized wave for Ed ^ 3. As Ed increases from unity, the phase speed Vph is practically constant for all 
s-polarized waves as well as in set pi; Vpu is almost constant in set p 2 for Ed ^ 1.5 and then increases linearly 
with Ed] and Vph is almost constant in set p 3 for Ed "^2. and then increases linearly with Ed- 

The propagation distances Ap^op of all three s-polarized compound SPP waves are almost equal when 
Ed = 1, and all decrease almost linearly with the same rate as Ed increases. All three p-polarized compound 
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Figure 4: Spatial variations of Px{0,z) (blue solid lines) and Pz{0,z) (red dashed lines) with respect to z 
for the compound s-polarized SPP waves, (a) L = 30 nm or (b) L = 200 nm in set Si; (c) L = 40 nm or (d) 
L = 200 nm in set S 2 ', and (e) L = 50 nm or (f) L = 200 nm in set S 3 . 
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-p-polarization 
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1.5 


2.5 
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Figure 5: Variation of Re(g) as a function of the relative permittivity Ed of the HID material when L = 25 nm. 



Figure 6 : Variation of the propagation distance Aprop 
material when L = 25 nm. (a) p-polarized compound 



as a function of the relative permittivity Ed of the HID 
SPP waves, (b) s-polarized compound SPP waves. 


SPP waves are different when Ed = I; whereas Aprop decreases almost linearly in set p 2 as Ed increases, the 
sets Pi and ps show more complicated dependencies of Aprop on Ed- 

Finally, Figs. andshow the spatial profiles of the Cartesian components of the time-averaged Poynting 
vector for two compound SPP waves belonging to the sets pa and S 3 , respectively, for two different values of 
Ed- Due to the significant coupling between the two dielectric/metal interfaces, a fraction of the energy of 
all compound SPP waves resides in the HID material and that p-polarized compound SPP waves are more 
sensitive than s-polarized compound SPP waves to that material’s relative permittivity. 


4 Concluding remarks 

We formulated and solved the boundary-value problem for compound surface waves guided by a thin metal 
slab sandwiched between a homogeneous isotropic dielectric material and a periodically multilayered isotropic 
dielectric material. Solutions were found to exist for compound surface-plasmon-polariton waves of both p- 
and s-polarization states at a fixed frequency (or free-space wavelength, equivalently). 
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Figure 7: Spatial variations of Px{0, z) (blue solid lines) and Pz(0, z) (red dashed lines) with respect to z for 
two compound p-polarized SPP waves in set ps when L — 25 nm and (a) = 1 or (b) Sd = 4. 




Figure 8 : Spatial variations of Px{0, z) (blue solid lines) and Pz(0, z) (red dashed lines) with respect to z for 
two compound s-polarized SPP waves in set S 3 when L = 25 nm and (a) = 1 or (b) £d = 3. 
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For any thickness of the metal layer, numerical investigations indicate that at least one compound surface 
wave must exist. It possesses the p-polarization state, is strongly bound to the metal/HID interface when the 
metal thickness is large but to both the metal/HID and metal/PMLID interfaces when the metal thickness 
is small. When the metal layer vanishes, this compound SPP wave transmutes into a Tamm wave. The 
behavior of this compound SPP wave clearly shows the coalescence of two surface-wave phenomenons that 
otherwise appear completely different. 

Additional compound SPP waves exist, depending on the thickness of the metal layer and the relative 
permittivity of the HID material. Some of these are p polarized, the others being s polarized. All of them 
differ in phase speed, attenuation rate, and field profile, even though all are excitable at the same frequency. 
We conjecture that the multiplicity is enhanced when the metal layer is thin if the relative permittivity 
of the HID material is smaller than the spatially averaged relative permittivity of the PMLID material. 
Although not explored here but extrapolation from previous work on SPP waves guided by a metal/PMLID 
interface [9 -11 13 , we can state that the period and the composition of the PMLID material also play a 
significant role in the multiplicity of compound SPP waves. The multiplicity and the dependence of the 
number of compound SPP waves on the relative permittivity of the HID material when the metal layer is 
thin could be useful for optical sensing applications. 
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